arXiv:1505.02060vl [quant-ph] 8 May 2015 


Time delayed control of excited state quantum phase transitions in the 

Lipkin-Meshkov-Glick model 

Wassilij Kopylov^ and Tobias Brandes^ 

^ Institut fiir Theoretische Physik, Technische Universitat Berlin, D-10623 Berlin, Germany 

(Dated: May 11, 2015) 

We investigate the role of dissipation in excited state quantum phase transitions (ESQPT) within 
the Lipkin-Meshkov-Glick model. Signatures of the ESQPT are directly visible in the complex 
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I. INTRODUCTION 

The properties of a quantum phase transition (QPT) 
[1] and an excited state QPT (ESQPT) [ 2 , 3] driven by 
quantum fluctuations in many body quantum systems at 
zero temperature feeds the interest since many decades. 
At a quantum level, the ESQPT is hidden in the sys¬ 
tem’s density spectrum as discontinuities or divergence 
[ 2 , 4]. For some famous models like Dicke superradiance 
[5] or Lipkin-Meshkov-Glick (LMG) [ 6 ], the main prop¬ 
erties of both QPTs can be explained and understood at 
a semiclassical level in the thermodynamic limit, where 
a QPT corresponds to a bifurcation [7, 8 ] and a ESQPT 
is connected to a saddle point in the semi-classical en¬ 
ergy potential [ 2 ]. Furthermore, the density of states can 
then be calculated analytically [4, 9]. For both models, 
the ESQPT manifests in the observable average as a peak 
at a certain energy [4, 10, 11] and was already observed in 
molecular systems [ 12 ] or microwave Dirac billiards [13]. 

Feedback control is a promising tool to change the sys¬ 
tem dynamics in a desired way. Influence of the laser 
statistics [14], neurosystems [15] or even control at a 
quantum level [16-20] are only some examples of its pow¬ 
erfulness. One usually distinguishes between the closed 
and open control loops, in the last one the feedback de¬ 
pends on the state of the system, where time delayed 
Pyragas control [ 21 ] has an important niche. For in¬ 
stance, it was used to speed up the convergence to a 
steady state in the dissipative quantum system [ 22 ]. 

In our previous works, we applied time delayed Pyra¬ 
gas control to a Dicke model [5] to create new non¬ 
equilibrium phases [23] and suggested a new method to 
extract the ESQPT-signal from the time evolution in the 
closed LMG system [ 6 , 24]. Furthermore, closed loop 
control was already applied to the LMG model to in¬ 
duce new phases or to modify the divergence in density of 
states [25-27]. Moreover, Oberthaler et al. implemented 
the LMG model using the existing ingenious experimen¬ 
tal setup based on Bose-Einstein condensates [29], 
which offer a high degree of freedom for system prepa¬ 
ration, providing a possibility to test new insights about 
the (ES)QTP. Surprisingly, the effects of a dissipative en¬ 
vironment on the ESQPT seems not to be well studied in 


contrast to the QPT [30], though they are always present 
in experimental realisations. 

Inspired by this, we study the effects of dissipation on 
the ESQPT signal in the LMG model and apply a time 
delayed feedback scheme to cancel them by the creation 
of new phases. In our model, we condition the atomic 
coupling on the difference of a spin observable average 
at two different times and perform the calculation at a 
mean-field-level. On top we show, that in dissipative 
systems the ESQPT is directly visible from a complex 
system spectrum. 

Our work is organised as follows. In section II we in¬ 
troduce the dissipative LMG model and the feedback. In 
section III we study the dissipative effects on a ESQPT at 
a quantum level evaluating the complex spectrum of an 
effective Hamiltonian. In section IV we show the smooth¬ 
ing effects of the ESQPT signal due to dissipation at a 
mean-field level and show in section V how to compen¬ 
sate them using time delayed Pyragas control. In the last 
VI section we discuss the results. 


II. LMG MODEL WITH DISSIPATION AND 
FEEDBACK 

The LMG model [ 6 ] describes an interaction between 
N spins and represents a special case of a Heisenberg 
Model. In general, the LMG Hamiltonian reads 

H = ( 1 ) 

where R = ^ collective angular 

momentum operators, h is an effective parameter for ex¬ 
ternal magnetic field in ^-direction and 7 ^, or 7 ^ describes 
the spin-spin interaction strength, which is the same for 
all spins in the LMG model. In the following, we will al¬ 
ways use the isotropic 7 ^ = 0 case, if we do not explicitly 
point to 7 ^ 7 ^ 0 case. 

Especially in the experimental realizations, the LMG sys¬ 
tem is always coupled to an environment which cause the 
damping and thermalisation of the system and can be 
modelled by a master equation [30, 31] with collective 
decay [32] 
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'p=-i[Hrp]-i^V[J^]p, (2) 

where V[J^]p = J~J^p + pJ~J^ — 2J+pJ“ is the 
Lindblad-Dissipator and = Jx ±iJy. 

To compensate the dissipative effects we assume a time 
dependent coupling jx of Pyragas form [21] . Therefore 
we condition jx to depend on a difference of averages 
at two different times with time delay r as following 

= = (3) 

In the following, we investigate the ESQPT signal at 
the quantum level using the effective Hamiltonian ap¬ 
proach and at a semiclassical level using the solution of 
mean-field equations in thermodynamic limit. In the sec¬ 
ond case we show the unique effect of the feedback loop 
especially in context of the ESQPT signal. 

III. DISSIPATIVE ESQPT AT QUANTUM 
LEVEL 

The fate of the ESQPT signal in presence of dissipation 
is still not studied well. In closed systems the ESQPT is 
hidden in the energy spectrum or is visible in observable 
averages as a function of energy [10] as a non-analyticity, 
which can be obtained by quantum mechanical or semi¬ 
classical calculations. Its origin can be understood at 
a semi-classical level as critical points from the energy 
surface. 

In dissipative systems described by a Lindblad master 
equation, an effective non-hermitian Hamiltonian can be 
defined in a standard way. Rewriting the master equation 
Eq. (2) as 

p=-i [Heff ,p] + -^ J^pj ~, (4) 

we obtain the effective non-hermitian Hamiltonian 

Heff = H-i^J-J+ (5) 

with complex spectrum, which is shown in the Eig. I 
(upper) for different values of k in the symmetry broken 
phase. Eor k 0 the spectrum has an imaginary part 
which scales with k. The imaginary part can be inter¬ 
preted as a decay rate at a certain energy level [33, 34]. 
To complete the spectral information we show in the 
lower part of Eig. I the corresponding density of states. 
Eor example, open triangles show the density of states 
along the blue dashed line in the upper figure in the 
known k = 0-case [9], where a logarithmic divergence 
at E = —0.5V is due to the ESQPT. 

How does the dissipation affect the ESQPT? The ES¬ 
QPT survives and, somewhat surprisingly, it becomes 
visible not only in the density of states of the non- 
hermitian Hamiltonian Vg//, but can be also directly 
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FIG. I. (upper) The complex spectrum for different values 
of 7 . The T-structure indicates the position of the ESQPT. 
The inset shows the zoom around the ESQPT energy, (lower) 
The density of states of the effective LMG Hamiltonian shows 
a logarithmic divergence at E = — 0.5iV . For comparison, 
the 7 y 7 ^ 0 case is additionally shown in both figures (open 
circles). Parameters: (7 = 1.5h),N = 1000 — 10000,A = 0. 


seen from its complex eigenvalues, see Fig. I (upper). 
The ESQPT is hidden in this representation in a Y-form 
at the energy of E = —0.5V. Another feature (due to 
the assumed Lindblad operator (2)) is the vanishing of 
the imaginary part of H^f f at the north and south poles 
of the corresponding Bloch sphere, which leads to zero 
imaginary part at the corresponding energies of ^0.5V. 
Thus the decay rate is zero there and the dissipative ef¬ 
fects disappear at this points. Note, this effect is also 
present at the level of the mean-field Eqs. (7), where dis- 
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sipative terms do not contribute at the poles for arbitrary 
n values due to conservation of the spin length. 

We used two different methods to calculate the density 
of states in the n = 0.05h case. First, we used only the 
real part of the eigenvalues and counted their number 
in a certain energetic window (blue crosses in Fig. 1). 
Second, we used 

u{E) = -I/tt^IY ( 6 ) 

with the non-hermitian Hamiltonian (open squares). We 
emphasize, that the results agrees and have still a loga¬ 
rithmic divergence at the energy E = —0.5N which can 
be well fitted by a log-function (red dotted curve). Espe¬ 
cially for E > —O.bN there is no strong deviation from 
the K = 0 case (open triangles). Only for E < —O.bN, 
there is a deviation from the non-dissipative result, which 
can be better seen for the curve with a bigger dissipation 
rate k = O.b (green dot-dashed). 

We emphasize, that only the Y-structure and not the 
zero imaginary part in the complex spectrum indicates 
the ESQPT. Eor 7 ^ = 0 both effects are at the energy 
E = —0.5, though they can be easily separated for an 
anisotropic LMG Hamiltonian ( 7 ^ 7 ^ 0), where the ES¬ 
QPT can be shifted to other energies and a jump in a 
density of states can occur on top [9]. Setting 7 ^ = 2.5h, 
we observe a shift of the Y-signal to the corresponding en¬ 
ergy of ESQPT in this case (open circles in the upper part 
of Eig. 1 ). The corresponding density of states (open 
circles in the lower part of the figure) has a peak at the 
ESQPT energy E ^ —0.54A^ and a jump at E = —0.5A^ 
[9]. Note, that the finite size effects smooths the peak 
and the jump in this case. 


IV. DISSIPATIVE ESQPT SIGNAL AT 
MEAN-FIELD LEVEL 

A. Mean-field-equations 

To derive a mean-field equation for the dissipative 
LMG system, we use = Tr{Op)^ assume the fac¬ 
torisation assumptions ( 0162 ^ ~ ^ 02 ^ which is 

known to hold in the thermodynamic limit and to forecast 
the same observable averages and the phase transition as 
the quantum mechanical calculations [24, 35]. We then 
obtain a following set of closed semiclassical equations of 
motion [30] 


jx{t) = hJy{t) - (7) 

jy{t) = -hJx{t) + - I^Jy{t)Jz{t), 

jt{t) = -‘^^x{t)Jx{t)Jy{t) + {Jx{tf + Jy{tf) , 

with rescaled averages Ji = Without time de¬ 

lay (A = 0) the QPT, which is one important property 


in this system, corresponds at this semi-classical level 
to a pitchfork bifurcation [31]: Eq.(7) has two station¬ 
ary states J?, corresponding to a normal phase with 
(J^, J^, J^) = ( 0 , 0 , 1 / 2 ) and a symmetry broken phase 


= ±1 


' —4h^(7 — ^ 7 ^ — + /^^(7 -h ^ 7 ^ — 


87/^^ 


( 8 ) 


J 


0 _ 7 - tO tO _ 7 - ^7^ - 


y 




whose stability swaps at a critical coupling 


lc-h + ^. (9) 

Note, that even the dissipative model still fulfils the con¬ 
servation law = 1/4, thus the dynamical evo¬ 

lution is restricted to a sphere. Furthermore, even with 
time delay (A > 0) the fixed points remain the same, as 
the Pyragas term vanishes in the steady state. 

In the following, we investigate the ESQPT signal in 
presence of damping with and without control. Later we 
show that the time delayed coupling 7 (t) may affect the 
linear stability of fixed points and the dynamical evolu¬ 
tion of the system in a completely unexpected way, par¬ 
ticular in its acting against the dissipation. 


B. Dissipative damping of the ESQPT 

Using the semiclassical equation of motion Eq. (7) we 
now study the action of dissipation on the ESQPT at this 
level. Therefore we first look at a dynamical evolution of 
the LMG system. The spin averages of the system are 
restricted to the Bloch sphere, which is shown in Eig. 2 
(inset). The color represents the energy for a given sys¬ 
tem configuration for /^ = 0 and the white lines represent 
the paths with the same energy in the symmetry bro¬ 
ken phase. Without dissipation the system follows one 
of this paths keeping the energy fixed, i.e. staying in the 
eigenstate. Eor each eigenstate one can compute an aver¬ 
aged value for (7^,7^). The black (dotted) curve (Eig. 2) 
shows this values for multiple eigenstates. This is a novel 
representation of two observable averages, which was re¬ 
cently suggested to visualize the energy independent ES¬ 
QPT signal [24], which is visible again as a peak. The 
continuation of the peak would end by (0,1/2). In the 
upper half of a Bloch sphere we see a separatrix, a (white) 
path which goes through a north pole of the sphere. Due 
to the symmetry of the path, its averaged values for 
and Jy should be zero, whereas the Jz average is 1 / 2 . 
Thus, paths close to the separatrix are responsible for 
the ESQPT peak in (7^,7^) diagram. Without dissipa¬ 
tion {k = 0 ) the period length for fixed energy (Eig. 2 , 
upper right) diverges at the separatrix energy [ 2 ]. 
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FIG. 2 . The figure shows the smoothing impact of damping 
on the ESQPT signal for different n values. The (blue) crosses 
show the averaging results using the integration method for 
At = 20. The ESQPT-signal under influence of damping gets 
smoothed. (Inset) Trajectories on the Bloch sphere visualize 
the system evolution Eq. (7) for Hi — O.l/i (dashed green), 
Hi = 0.2/i (orange) and = 0 (white contours). The color 
shows the rescaled energy in /c = 0 case. Parameters: 7 = 
1.5/1, A = 0. 


For Hi the energy is not conserved any more and the 
system tends oscillating around the J; 2 -axis to a steady 
state Eq. ( 8 ). Fig. 2 (inset) shows two examples of the 
system state evolution {Jx^ Jy^ Jz){t) for fwo different k 
values, which where obtained by solving Eq. (7) numeri¬ 
cally. The ESQPT signal is now hidden in the dynamical 
evolution of the system. But, as we can see in Fig. 2 , 
especially for big k values there are only less paths, which 
are close to the separatrix, thus the ESQPT signal will 


A = -4A J° 


( 4 ) 


jO jO 
^ X ^ y 






^ Jz 

•27^ 


with eliminated Jz component by use of the spin-length 
conservation law. 

The roots of the corresponding characteristic equation 

det(Al — B — Aexp(—Ar)) = 0 (13) 

determines the stability of a fixed point, which is stable 
if all real parts of all solutions A are negative [37]. 

In Fig. 3 we plot the biggest real part of eigenvalues 
in the r — A— domain. We see that there is a window 


be damped especially for big dissipation rates k. The 
impact of damping to the ESQPT signal in the 
plane can be obtained by calculating the mean values of 
Jj(t), Jz{t)-eYo\ution for a period or by finding some op¬ 
timal effective period At, in a way that the mean values 
calculated using the definition 

0 {t) = ^ 0 (t)dt, a G { 0 , Tmax}. ( 10 ) 

matches to the closed case as good as possible. Note, 
as a period length in the first case we take the time for 
one full rotation around the T^-axis, which changes with 
time. 

In the first case, the results for different damping rates 
hv are shown in Fig. 2 by colored circles/diamonds and 
unfilled squares. We see, that the peak is now smoothed, 
but still visible especially for small damping rates. The 
blue crosses show the second case with effective period 
which lead to a much better ESQPT signal. 

V. PYRAGAS CONTROL OF THE ESQPT 
SIGNAL 

Now we set A 7 ^ 0 in jx{t) Eq. (3) and investigate the 
impact of time delayed control on the system. 


A. Linear stability analysis with time delay 

A usual approach to analyse the effects of time delayed 
feedback is first to check the stability of fixed points in 
the presence of control [36]. Therefore we linearise Eq. 
(7) around the fixed points Ji{t) = Jf + 6 Ji, obtaining 
the following system of linearised equations with Sv = 


dtSv{t) = B • 6 v{t) + A • 6 v{t — r), (11) 

with 


tO tO \ 

- kJO h + \ 

+ -h + 27 JO kA- + 4AJ°7 - 27^^ - k4 j ’ 

( 12 ) 


for A -values, there the stability of fixed points oscillates 
from stable to unstable and from unstable to stable while 
increasing the time delay. Outside this window, the fixed 
points remains either stable (A ^ Ih) or lose their sta¬ 
bility forever (A > 2.5h). Note, that the boundaries be¬ 
tween stable and unstable zones (brown line) can be cal¬ 
culated from an analytical expression, see Appendix A. 
Next, we analyse the system properties in the unstable 
regime, use them to obtain a sharp ESQPT signal and 
show chaotic behaviour for larger time delays r. 
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FIG. 3. Stability diagram in the r — A-plane. The color 
represents the stability robustness of fixed points in the sym¬ 
metry broken phase, which are not stable for positive val¬ 
ues (red coloured area). The brown line defines a bound¬ 
ary condition and has an analytical expression. Parameters: 
7 = 1.5/i, K = 0.05/i. 


B. Feedback compensates dissipation 

Our feedback scheme Eq. (3) modifies the system dy¬ 
namics in an interesting and unexpected way. Increasing 
the time delay r, we cross the boundary and make the 
fixed point unstable. For smaller r values the trajectory 
ends in a new stable state in form of a limit cycle, thus a 
Hopf bifurcation occurs. Moreover the trajectory of the 
limit cycle has only small deviations from paths with 
fixed energy, which the LMG system would take without 
dissipation. The size of trajectories can be changed by r, 
thus tuning the time delay value corresponds to a change 
of energy in a closed LMG system. Fig. 4 demonstrates 
the feedback action, showing the trajectory evolution 
for different values of r. Here we start close to the fixed 
point. The red (thick) curve shows the stationary state. 
Note, that the change of initial condition to the region 
outside the separatrix leads to the same stationary state 
for the considered r values. 

Thus, using the stationary limit cycle states for fixed 
r values offers a possibility to obtain an ESQPT-Signal 
again, by calculating the 

j,=-^ r Jz{t)dt ( 14 ) 

^2 “ Jti 

and (with similar definition) averages for fixed val¬ 
ues of r (t 2 > ti and with ti big enough to become a 
stationary solution). Fig. 5 shows the results of this cal¬ 
culation and compares them with the closed and dissipa¬ 


tive cases without feedback. The (orange) rhombi shows 
the Jz , Jx time averages in stationary limit cycle phase 
of the dissipative LMG system, the black curve shows 
the ESQPT -signal of a closed system without control. 
Each rhombus has its own time delay, the inset in Fig. 
5 shows the r-dependence of the averaged values. Note, 
that up to r ~ 0.2/h the fixed point is stable and the 
averaged values are the values of the fixed point. We see 
a very good overlap between the control caused and the 
original signals, thus our feedback schema compensates 
the dissipative smoothing (unfilled green squares) of the 
ESQPT-signal very well. 


C. Chaotic behaviour 

For r ^ 1 the stationary dynamics can become much 
more complex than just a creating limit cycles. Oscilla¬ 
tions with more than one maximum and minimum ap¬ 
pear, which is known as a way to chaos by period dou¬ 
bling [38, 39]. In Fig. 6 (upper) we plot the maxima and 
minima of Jx{t) oscillations for t ^ 1 as a function of 
time delay r and mark different areas by capital letters 
A-E. In the inset we show a zoom for the area A. Increas¬ 
ing r from zero, the stable fixed points (arrow) loses their 
stability (orange dotted line) at the first boundary condi¬ 
tion from the Fig. 3 and a Hopf-bifurcation appears (part 
A). In part B, the fixed points become stable again. But 
there exist still a stable limit cycle solution with rather 
big Ja,-amplitude. Thus, in area B two stable solutions 
are possible, which is not in contradiction to the stability 
of the fixed point, as the initial condition is not chosen in 
a way to fulfil the linearised assumption. For 4 < r < 5 
this limit cycle disappears. Further increase of time delay 
leads to a creation of a period doubling structure (area 
C), which is separated by windows where the solution 
converges to a limit cycle. In the region D the fixed point 
becomes stable again and the double period structure is 
gone, whereas in the Region E it appears again. Note, 
that the dotted line at = 0 represents the unstable 
trivial fixed point. 

Such chaotic behaviour can also be used to obtain the 
ESQPT signal, see blue crosses in Fig. 5. Fixing the time 
delay in the chaotic area of the phase diagram (Fig. 6), 
one can use the integration method with effective period 
Eq. (10) to obtain the shown curve from the dynamical 
system evolution on the Bloch sphere. We see, that also 
this result matches pretty well to the original ESQPT 
signal. 


VI. DISCUSSION 

In this paper we have demonstrated the effect of dis¬ 
sipation on the ESQPT signal for the LMG model and 
showed how to compensate it using time delayed Pyragas 
feedback modulating the interaction parameter between 
the atoms. Our results show, that the ESQPT is encoded 
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FIG. 4. Dynamical evolution of the controlled LMG system on the Bloch sphere for the following r-values (from left to right): 
rh = 0.2; 0.25; 0.3; 0.31; 0.35; 0.5; 1; 2.5; The red thick path represents the stationary state. Increase of r forces the solution to 
cross the separatrix. In this way the ESQPT signal (Fig. 5) is restored from stationary solutions with different r values. Note, 
the axes and scaling are as in the inset of Fig. 2 but are not visible due to breakdown. Parameters: 7 = 1.5/i, A = 1, = 0.05/i. 



FIG. 5. Time delayed feedback undo the dissipative effects 
to the ESQPT signal. The orange (light-coloured)rhombus 
shows averaged system values Jz , Jx in the stationary state 
for hxed r with r G [ 0 . 1 , 0 . 6 ]/i, which matches to the sig¬ 
nal of the closed system (black curve) very well. Unfilled 
squares represent the ESQPT-signal under dissipation. The 
blue crosses shows the ESQPT-signal obtain from a dynami¬ 
cal evolution of the system for one big r-value using Eq. (14). 
(Inset) Jz, Jx as a function of r in a steady state. Parameters: 
7 = 1.5/i, A = Ih. 


in the spectral properties of the effective Hamiltonian as 
well as it is visible in the measured averaged values of 
the spin components. In the last case, smoothing effects 
appear which can be undone by our feedback scheme. 

We think, that an experimental measurement of an 
ESQPT signal using the time delayed method is easier 
than in an ideal, closed system. Using time delay one 
has only to measure the system values for different time 
delays r, instead of preparing the system in eigenstates or 
coherent superpositions to obtain the same information 

[24]. 

We also checked the interaction strength 7 ^, depending 
on other operator differences instead of J^, or a modu¬ 
lated magnetic field h instead of 7 ^,. However, in both 
cases the general dynamical properties remains the same. 
The LMG system has still limit cycles close to the separa¬ 
trix for some fixed time delayed values and has a param¬ 


A B C D E 



EIG. 6 . Bifurcation and period doubling for different time 
delays r. Eor each r-values the values of local maxima and 
minima of Jcc(t)-value are plotted in stationary case. The 
dashed orange lines shows the unstable fixed points. The ar¬ 
rows points to the non-trivial fixed points. The phase diagram 
is divided the regions A-E with different properties. The inset 
shows shows a zoom of area A. Parameters: 7 = 1.5/i, A = Ih. 


eter range with chaotic behaviour. Both effects might be 
interesting from an experimental point of view. On the 
one hand, it is easier to control the magnetic field /i, on 
the other hand, choosing another feedback loop can shift 
the appearing effects to other r values, which could be 
easer to realize. 

We also tried to find ESQPT signatures in the corre¬ 
lation function [40] Ct{z) = ^ J+(t + using the 

quantum regression theorem [41], but have not succeed, 
as the ESQPT signal (which should be visible as a peak at 
zero frequency in the Fourier space) interfere with effects 
like macroscopic occupation and degenerated spectrum in 
the symmetry broken phase. However one possible way 
out could be to drive the LMG system additionally with 
an external laser and calculating the resonance fluores¬ 
cence spectrum, but this would require a re-derivation of 
the master equation. Note, that such spectra have been 
already calculated for an LMG system, but only in the 
linearised version [31]. 

Similar to the Dicke model with feedback [23], the 
Pyragas controlled LMG model has stable limit cycle 
phases. In contrast to the Dicke model, it shows chaotic 
dynamics for bigger time delay. This is surprisingly as 
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one would expect such behaviour especially from the 
Dicke system, which is chaotic in nature [42]. How¬ 
ever,the Dicke model has an additional bosonic mode 
which is not bounded by a conservation law like the spin 
components. This could be a reason, why no chaotic be¬ 
haviour appears there. 

We think, that the feedback-induced limit cycles are 
hidden property of the LMG modell. On the one hand 
limit cycles are a natural property of the closed system. 
On the other hand, the time dependent limit cycles de¬ 
scribe the ESQPT signal pretty well. Though the chaotic 
behaviour is an artefact of the time delay feedback, as it 
was neither a part of a closed LMG system. 

The feedback scheme is applied at a semiclassical level, 
thus we have neglected the influence of fluctuations in 
the thermodynamic limit, which could be important [43]. 
Nevertheless, we think, that the fluctuations have not 
dramatic contributions to the described effects for N ^ 1 
as they should scale with 1/\/N and the used feedback 
scheme does not modify 7 in a strong way. Furthermore, 
the semiclassical LMG model predicts in many cases the 
same results [4, 10, 24]. But the full quantum version of 
the considered feedback type still remains an open issue. 
However, a recently published article [20] shows a way 
to go beyond mean field for a coherent type of feedback 
where the author describes feedback action via mapping 
to a bigger system. This would be one possible way to 
study the role of oscillations in quantum systems with 
one special feedback type. 
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Appendix A: Boundary condition 

To determine the boundaries in Fig. 3 we choose only 
imaginary A values in Eq. (13), thus A = i • 5,5 G [R. 
Splitting then the Eq. (13) in an imaginary and real 
part, we obtain the correspondingly conditions 


G , = h ^- S ^ + - 27/17 - 

+ - 4/i(7)V7 - 4^77(7)7, 

G 4 = 4/i(7)"7a + 4^77(7)7, 

G5 = 4777AS. (A3) 

Bringing all sin and cos terms in both equations to one 
side, squaring them and adding together, we eliminate 
the r-dependence and obtain the following equation 

0 = + Fis^ + Fo, (A4) 

with 

^’o = (/iV° + 27 /i (( 7 )' - ( 7 )') 

+ «27(-(7)^-(7)V(7)^)^. 

(^/i 27 + 2h (( 7 )' (7 - 4(7)7) - 7(7)') 

+ «7 (-77)' - 877(7°)'a + AC ((7)' - (7)')) 

)}■ 

= 7^ {“^^'(7°)' + 4/1 ((7)' (2(7)'a - 77) + 7(7)') 
+ - 4«(7)7° (7 - 2(7)7) 

+ 2{4f [{J^f + 27 (7 - 4(7)' a )) - ac 2(7)^) 

- 4^77 ((7)' (7 - 2(7)7) + 277 - 27(7)') 

+ ((7)" - 2(7)^(7)^ + 2(7)") |, 

which can be solved for s and has then in general 4 dif¬ 
ferent solutions 


0 = Go + Gi • cos { st ) + G 2 • sin(5T), 
0 = G3 + G4 • cos{sr) + G5 • sin(5T), 


Go = - 


77)'a 


7 

«(7)'5 


2jJ^J^s non 

-4777Ag 


“1“ 2kJ2s, 


7 

Gi = 4777 A 5 , 

G 2 = -4/1(7)77 - 4^77(7)7, 


(Al) 


(A2) 



The Eq. (A5) fixes now the eigenvalue A = is for 
a given fixed point and feedback strength A. Eq. (Al) 
gives for every fixed 5 -value the corresponding time-delay 
r, solving for example the imaginary part for r we obtain 


1 

T = - arctan 
5 


g 2 + g 2 ^g 2 +g 2 ^ 

Gi 


iG 0 G1 
Gf+Gl 


+ 


27r 

5 


(A6) 



where z ^T. The choice of z is necessary to get bound¬ 
ary conditions at higher time delays r. Note, that the s 
dependence is also hidden in Gi. 

At this step we have still to much solutions. A lot of 


them are non-physical (if r < 0 or r G C) or do not fulfil 
the real part equation Eq. (Al) and have to be sorted 
out. The remaining solutions are plotted as a brown line 
in Fig. 3. 
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